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Functional relations play a key role in the study of integrable models. We argue 
in this paper that for massless field theories at zero temperature, these relations can in 
fact be interpreted as monodromy relations. Combined with a recently discovered duality, 
this gives a way to bypass the Bethe ansatz, and compute directly physical quantities as 
solutions of a linear differential equation, or as integrals over a hyperelliptic curve. We 
illustrate these ideas in details in the case of the c — 1 theory, and the associated boundary 
sine-Gordon model. 
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Integrable models are one of the most widely studied areas of mathematical physics. 
The fact that there are an infinite number of conserved charges commuting with the Hamil- 
tonian means that in principle physical quantities can be computed exactly. In practice, 
this task is usually quite difficult. 

Indeed, there are surprisingly few techniques at our disposal. One of the most suc- 
cessful ones is actually a related set of techniques going by the name of the Bethe ansatz. 
In situations where the Bethe ansatz is applicable, one makes some technical assumptions 
about the wave function, and then can derive a set of equations for the zeroes of the wave- 
function. The number of such zeroes goes to infinity in the continuum limit, and one can 
then use various density functions and the integral equations they satisfy. The equations 
obtained are rarely solvable in closed form but a great deal of quantitative information can 
be gained by studying them in various limits or by numerical solution. 

Since the Bethe ansatz is usually technically formidable and often intractable, it is 
quite desirable to develop an alternative formulation. One method being explored in recent 
years is the use of functional relations or "fusion" relations for physical quantities. The 
existence of such relations is almost obvious, and based on very general properties of 
the Yang-Baxter equation and quantum groups. The assumptions made generally involve 
various analyticity properties of physical quantities, and can often be verified using the 
other approaches. Quite frequently, it can be shown that the Bethe ansatz equations follow 
from the fusion relations, but the latter seem more powerful, more convenient, or deeper 
i- 

Another method that appears promising is the use of duality. It has long been known 
that some quantities in certain integrable boundary field theories (for instance, the current 
in the boundary sine-Gordon model [|2|, as well as the magnetizations in the Kondo model 
0) do satisfy a weak to strong coupling duality. Moreover, inspired by work in supersym- 
metric gauge theories, it was shown in |@J|,H that these quantities can be expressed, at 
T = 0, in terms of integrals on hyperelliptic curves. It would be quite useful to understand 
the precise mathematical framework underlying these properties. It seems possible that 
when combined with some sort of analyticity and maybe some other type of information, 
they can be used to bypass the Bethe ansatz, and lead directly to the determination of 
physical quantities in a way analogous to the supersymmetric gauge theories, as pioneered 
in0. 

In this paper, we present some further progress along these lines, restricting to the 
simplest case of the boundary sine-Gordon model, and the underlying c = 1 conformal 
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field theory. At zero temperature, we show that it is possible to interpret the fusion 
relations as monodromy relations for the current (and magnetizations). This, combined 
with duality and an assumption about the number of singularities in the plane of complex 
couplings, leads to a direct determination of physical quantities that indeed bypasses the 
Bethe ansatz. 

A key ingredient in the following is the current. Although it was introduced on purely 
physical grounds in , this quantity turns out indeed to also be central to our alternative, 
quite mathematical, approach: in a nutshell, we will see that the current is the unique 
solution of a differential equation that is regular at the origin, while for instance densities 
of physical excitations can be expressed as monodromies around the other singularity (1 or 
oo). The current is so special because it is exactly self-dual. Intuitively, the property arises 
because the current is measuring charge transport, so only operators with charge affect its 
computation (in contrast with quantities like the free energy or the magnetizations), and 
they are strongly constrained by integrability f^. 

For convenience, we summarize here the notations to be used in the following: 



with A defined so that z{9 = A) = —1. 

1. Analytic properties of the current 

Since in the case discussed in this paper, we do have the Bethe ansatz solution avail- 
able, it is easiest and most concrete to use it to first illustrate the monodromy properties 
and the result. In the next section, we will show how all of these results follow without 
recourse to the Bethe ansatz. 

1.1. The model and the current 

The field theory discussed in this paper is usually known as the boundary sine-Gordon 
model in its massless limit. It consists of a 1-dimensional boson (p{x, r) on the half line 
a; > 0: there are no interactions in the bulk and the action is the usual (t/An) J dxdr^d^cf))'^ , 
but on the boundary there is an interaction Acos[(^(0)]. In addition, there is a voltage 
V coupled to the U{1) charge. The dimensionless couplings in this problem are u oc 
(^tiie coefficient of proportionality is non-universal and of little interest here) 




u = 
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and t (in the notation of m = e^+^-^B g^j^^ f _ l/gy Although physical quantities are 
continuous in t, many of the formal properties of the system are much easier to understand 
at rational or integer values of t. This is familiar from the Bethe ansatz study of the XXZ 
spin chain, which, with appropriate boundary conditions, is a lattice discretization of this 
model. 

The interaction violates charge conservation, so one can define the "current" describ- 
ing this charge violation. This situation is discussed at length in 0, where integral Bethe 
ansatz equations are derived for this current. Here we work with the normalized quan- 
tity X(u) = tl/V. At T = 0, the integral equations become linear, and the Wiener- Hopf 
technique allows one to find an integral expression for the current. This integral expres- 
sion yields explicit weak-coupling and strong-coupling series expansions. In terms of the 
coupling the expansion of the current for u large, real and positive (the UV) is 
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J(^.) = l-}_^a.(^)^.-2"(*-^)/^ (1.1) 



and for u small (the IR) is 

oo 

Iiu) = J2^n{tW'^'-'\ (1.2) 
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where 



n 



2 r[f + n(t-l 



As noted in 0, a duality t 1/t is already apparent in the above expressions. It was 
shown in that these expressions for the current can be reformulated in terms of an 
integral on the hyperelliptic curve y"^ = x + — u^: 

l{u) = 1-^1- (1.4) 

where the contour C starts at the origin, loops around the branch point on the real axis 
(when u is real), and returns to the origin. This result is continuous in t, but already 
we see a difference between t rational or irrational: the curve is of finite genus when t is 
rational and infinite genus when it is not. This expression means that the current obeys a 
differential (Picard-Fuchs) equation, whose order depends on the genus of the curve. 



1.2. Generalized hypergeometric functions 

To discuss the extension of this current to complex couphngs, and the associated 
analytical properties, it is convenient to consider further the expansion ( |1.2| ) . For simplicity, 
we restrict to the simplest case when t is an integer. We set 



-e-2A(t-i)^2(t-i). 



(1.5) 



where 



-1 



[tlnt-(t-l)ln(t-l)]. 
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The expansion ( |1.2|) defines a unique analytic function for \z\ < 1. Introducing the gener- 
alized hypergeometric function 
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and using the duplication formulas for F functions, one finds the simple expression 
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For t = 2 for instance, we have 
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Like for usual hypergeometric functions, the analytic structure of the current in terms of 
the 2;-variable is quite simple. It has two singular points at ^ = 1 and z = oo , and there is 
a cut running from 1 to oo. The monodromies around 1 and infinity are identical, and the 
Riemann surface of the function has t sheets (except for t = 2 where there are logarithmic 
terms, as is obvious from ( p..8| )). 

We will write the current for \z\ > 1 as 



1 = 1 



n=l 



U 



-2n{t-l)/t 
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(1.9) 



4 



where < arg z < 2n. This expansion, which we obtained originaUy from the Bethe 
ansatz and then from duahty considerations, can also be deduced from standard references 
on the theory of generahzed hypergeometric functions (Meijer's G functions): see for 
more details. As an example, let us recall the well-known result for t = 2: 



1 - 14 ^) = l-r(3/2)r(l/2) {e-^-z)-'/' F (0, 1; I z-'^ 

r(3/2)r(-i/2 
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where in fact F (O, ^; ^; z~^) = 1. For t = 3, elementary calculations give: 
1 2 



;i.io) 



1 — 3F2 



3' 3' 

3 5 

4 ' 4 



1-ci (e-'"^)'/^ 



7 12 
12' 12' 3' 



-C2 (e-^"^)'^V ( —^—;Z.z- 
2 V J \ 12' 3' 



15 4 



C3 (e-'^'z) 3F2 



5 3 
4'4' 
5 4 
3 ' 3 



; z 



(1.11) 



where the numerical constants are 

r(i/3)r(3/4)r(5/4) 

r(5/i2)r(ii/i2) 
r(-i/3)r(3/4)r(5/4) 
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r(i/i2)r(7/i2) 
r(-i/3)r(-2/3)r(3/4)r(5/4) 
r(-i/4)r(i/4)r(i/3)r(2/3) ' 
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1.3. Differential equations 

Using standard results from the theory of generalized hypergeometric functions, it 
follows that for t integer, the current solves the following linear homogeneous differential 
equation of order t 



d f d 3 + 2k 

k=0 ^ ^ ' 



n(4+7 



(1-X)=0. 



(1.13) 



This can be verified using either (|1 . 7| ) or (|1.9| ). For example, for t = 3: 
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However, the expression ( |1.4| ) of X in terms of a hyperelliptic curve actually requires 
that X obey a differential equation of order t — \ (the Picard-Fuchs equation). This is easy 
to see: taking h derivatives with respect to u gives integrals of the form 



dx 



2b+l 



y 

where a is some integer. If a > t — 1 or greater, this can be reduced to a sum of integrals 
with a < t — 1 by integrating by parts. Explicitly, 



X 



c y 
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dx 
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X 



a—t 
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X 



a-t+l 



,26+1 



t(6- 1/2)^26-1 ' t{b- 1/2) 

The constant term arises because the contour is not closed: the integration by parts results 
in a non- vanishing surface term when a = t — 1. By using this relation, one can express 
the order t — 1 derivative of T as a linear combination (with u-dependent coefficients) of 
the lower derivatives, so therefore X satisfies a differential equation of order t — 1. Because 
of the constant term, this differential equation can have a inhomogeneous piece. 

Therefore ( |1.13|) must be a total derivative. This can easily be checked to be true 



directly: the reason can be traced back to the fact that the second product in (|1.13|) runs 
up to /c = t. Let us illustrate this for t = 3: ( |1.14| ) reduces to 



8z 1 , , . d , 
9 4 ^ ^ dz ^ ^ dz^ 



:i-X)=A 
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(1.15) 



The value of the constant At is determined by the normalization of X (here As = 1/4). 
The situation is similar for other values of t: the current is a particular solution of an 
inhomogeneous differential equation of order t — 1 which reads 



t-2 



.k=0 



3 + 2A; 
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(1-X)=A 
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(1.16) 
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k = l. 



t-l. 



One obvious question is the significance of the other t — 1 solutions of the differential 
equation ( |1.13|) . These solutions are also necessarily solutions of ( |1.16|) with At = 0. We 
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will see in the next section that, at t integer, they coincide with the densities of states of 
the various species of quasiparticles. 

Mathematically, it is convenient to express these other t — 1 solutions in terms of the 
functions 



Xk{z) = I 



J 



(1.17) 



defined for k = 1 . . .t — 1 and \z\ > 1. In other words, we take the difference of the current 
on two successive sheets. The functions Ik are defined for \z\ < 1 by first using ( |1.17|) for 
1^1 > 1, and then continuing around the singularity. The series expansions for z large are 



r(T) 
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\ r(n) r[| + n(i-l)] 
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Here it is assumed that < arg z < 2n. These coefficients differ from those in the 
expansion of the current only by the piece siii{nn {t — k) / 1) . Since the solutions of the 
differential equation for \z\ large are of the form 

z''/\ak{l/t) + ak+t{l/t)z + ...) 



it follows immediately that a basis of solutions of the differential equation consists of the 
t quantities I{z),It-2p{z) and Jt_2p-i {e'^'^z) = l[_2p_i{z). 

For the monodromies to be discussed in the next section, we will also need the IR 
expansion valid at small z. This does not follow instantly from ( |1.17|) , because the relation 
( p..l7|) does not hold at small z. These expansions must be determined by continuing the 
expansions in ( |1.18| ) around the singularity at 2; = 1 (this is easiest to do by writing a 
contour integral whose residues give the terms in ( p..l8| ); this form arises naturally in the 
Wiener-Hopf analysis). One finds 



int 



2{t 



-rruk ( 
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(2n+l)" 
L 2(t-l)J 
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{2n+l)t~ 
[ 2(t-l) J 



n=0 

2n + l 



n + 1 



cos 
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(1.19) 



g Z7rg2(t-1)A^\ 2(t-i) 



/C = 1, . 



where nik = 2 sin(/c7r/(2t — 2)), G+{0) and G+{i) are some numerical constants (see for 
their exact values). The function It-i is given by (|1.19| ) for k = t — 1 with an additional 
factor of 1/2 multiplying the right-hand-side. The radius of convergence of these expansions 
is the same as for the current: they converge for |z| > 1 and 1^1 < 1 respectively. 
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1.4- t rational 

The situation is rather similar for rational t, and depends mostly upon the value of 
the numerator Q int = Q/P. Here, we shall simply give the differential equation satisfied 
by the current: 



(1-J)=0, (1.20) 



where now 



- (1.21) 

This is a homogeneous differential equation of order Q. The fact that the current can be 
expressed in terms of the curve = -\-x^ — u^ means that this equation can be reduced 
to an inhomogeneous equation of order Q — 1, just as in the case t integer. 

1.5. Monodromies 



The differential equation ( |1.13| ) is analytic in z, so when one continues a solution 



around a singularity, one must obtain a solution upon return to the starting point. How- 
ever, the continuation of a given solution is not necessarily the same solution: it only 
need be some linear combination of all the solutions. This linear combination is called 
a monodromy. Studying the monodromies is useful for understanding the solutions; in 
fact, if one knows all the monodromies, then one can reconstruct the differential equation 
automatically. In this section, we construct explicitly the monodromies of the solutions 
from the small- and large-^ expansions. Later, we will reverse the order of the logic: we 
will derive the monodromies from general arguments and therefore infer the solutions. 

The monodromies around infinity give the behavior of the solutions when we move 
z — s> e'^'^'^z for fixed \z\ > 1. Since there are t solutions for the differential equation, 
the Riemann surface for the current is t sheeted, with a cut running from 1 to oo. The 
monodromies around infinity for the the whole set of solutions of the differential equation 
can be easily written by using the UV expansion. We have 

It-2p [e^'''^z) =Xt_2p+i -Tt-2p +I't-2p+i (1-22) 

^t-2p-l (e'^'^'z) =Xt_2p+i -Xi_2p -|-X^_2p_i -Xt_2p-2 +X^_2p_3 
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The last "boundary" terms depend on whether t is odd or even. For instance, when t = 2 

I (e^'^'z) =I(z) - 1\ iz) 

(1.23) 

X[{^^-z)^-X[{z) 

and when t = 3, 

l{^^-z) =X{z)-l'^{z) 

X'^{e^'^z)=-X^{z) (1.24) 
I^{e^'^z)=X'^{z)-Xt{z) 
The monodromies around the origin are also quickly worked out from the IR expan- 
sions. The current in particular is regular at 2; = 0, and we have 

X{e^'^z) =X 
T;_i (e2-z)=Xi_i-T,_2+X^_3 
Xt_2 {e^'^'z) =2X;_i +X;_3 -Xt_2 

1t-2p [e^^'^'z] =Xi_2p+i -Xt-2p +X[_2p_i 
1t-2p-l {e^'^'^z) =X^_2p^i - Xt_2p -\- l't-2p-l -1t-2p-2 +X^_2p_^ 



(1.25) 



(1.26) 



The last boundary terms depend on whether t is odd or even. For instance, when t = 2 

X {e'^'^'z) =X 
X[{e^'^z)=-X[ 
and when t = 3, 

I [f'^z] =X 

Xi {e^''' z) =2X'2 - Xi (1.27) 
X'2 (e^-^) =X'2 - Ji 

Finally, the monodromies around z = 1 follow from the other two, because going 
around a circle at large z is equivalent to going around z = and z = 1. Since the 
monodromies take a solution to a linear combination of the other solutions, we can repre- 
sent tlicm conveniently by matrices M^o, Mq and Mi. The conventions for orienting the 
monodromies are usually taken so that one has M^q = MqMi. For instance, one has for 
t = 3, 

/I 0-l\ 
Ml = Oil 

\o 1/ 
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One checks on this example the general property that the monodromies of the current 
around z = 1 and z = oo are identical, since the current has a trivial monodromy around 
the origin. 

We observe that one has to be careful in performing the continuation of the X[_2p_i 
quantities for \z\ < 1, since the expressions ( [1.18|) ,( |lT9| ) hold only for < arg z < 2%. 
But this continuation is easily worked out by using the previous relations. 

A striking feature of the relations (|1.22|) , (|1.25|) is that the monodromies for the den- 
sities do not involve the current itself. This is because (or conversely, requires that) they 
satisfy a homogeneous linear differential equation (|1.16| ) of order t — 1, one order less than 
the hypergeometric equation ( |1.13| ) we started with. This order t — 1 differential equation 
is the same one that the current satisfies, except that it does not have an inhomogeneous 
term. 



2. A determination of the current from duality 

In this section we will show that the monodromies derived above from the explicit 
Bethe ansatz solution can be obtained directly from general arguments, at least when t 
is rational. Therefore, the explicit solutions can be found without recourse to the Bethe 
ansatz. 

As was discussed at length in [p],[TO[|, the partition functions Zj of the spin j /2 Kondo 
model satisfy, as a consequence of the Yang Baxter equation, some general fusion proper- 
ties. There are various ways to write those. A particularly convenient one here is based 
on the quantities Yk = Zk+iZk-i; one has then, in standard notations 

Yk (g^/^x) Yk (q-i/^x) = [1 + lfc+i(x)] [1 -f- Yk-i{x)] , (2.1) 

where q = e*'^/*, and x oc A oc (which is oc e*^^~^)^^ in the notation of 0). For t 

an integer, these relations hold only for /c = 1, . . . , t — 3. Their closure is a slightly more 
delicate matter, that was discussed in detail in section 4 of [|T^. It requires the introduction 
of a pair of objects, Y±, for which one has = e^, and 

Yt-2 (?'/'x) Yt-2 (q-^/'x) = [l-f Yt-s{x)] [1 -t- y+(x)] [1 + y_(x)] , (2.2) 
together with 

Y+ (gi/^x) Y+ (^q-^/^x^ y_ (q^/^x) y_ (g-i/^x) = [1 -f Yt-2{x)f . (2.3) 
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Finally, we need the relation between the current and partition functions conjectured in 



(,--.)^-!^4.„i±^. (.4) 

This is true at any temperature T. 

The claim we are making is that the relations (|2.1| ) , ( p.2| ) , ( p.3|) together with (|2.4| ) , and 
duality, are equivalent to the monodromy relations studied previously, and hence lead to 
a complete determination of the physical quantities in the problem of interest. 

General considerations about perturbation near the UV and IR fixed point show that 
the quantities Y should have a singularity at x = and at a; = oo. In addition, there 
should be at least one other singularity. We have to make the hypothesis there is only 
one such singularity, which we set at z = 1 by convention. There are thus two cuts in 
the complex z plane for our quantities, one from — oo to the origin, and one from 1 to oo. 
The fusion relations are now simply interpreted as monodromy relations. To show this, 
we first observe that we can set Yk = exp (e/c is necessarily positive to guarantee that 
the partition functions are not trivial, even in that limit). Writing the relations in terms 
of our usual variable z, and continuing them to the entire complex plane, leads to 

ekie'^'z) + ek{e-'^z) = ek-i{z) + e^+M = 1, . . . , t - 3 (2.5) 

where we set cq = 0. Thus if oc -z^, these relations are equivalent to all but the first 
two monodromy relations in ( |1.22|) , ( |1.25|) . 

Next, we set Yj^ = e'^'"/^ and Y- = e*^"/^. It is not clear now what the sign of these 
e's is, since the Y± have no direct physical meaning as partition functions. To solve this 
difficulty, we need duality: the latter implies that the current should expand in powers of 
z'^ in the IR i.e. have no non-trivial monodromy around the origin. This implies that 
the right hand side of ( pj.4| ) vanishes, and thus, setting Y^ = e*^^/^ and YL = e^°-l'^ , that 
both eg and have a positive real part and are equal for 1^1 < 1. As a result, the fusion 
relations for < 1 close with 

Jt_2(e^"z)+Tt_2(e-*"^) =Tt_3 + 2Jt_i, (2.6) 

where 

1 / dts dta \ 

^ Recall that, as discussed in Q, this follows because the current is entirely determined by the 
leading irrelevant operator near the IR fixed point, unlike the magnetizations or l^'s. 
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and the coefficient of proportionality is tfie same as for X^. No special relation is required 
for the current this time, since it has no monodromy around the origin, 

T(e^^^^) = J. (2.7) 

Finally, the last relation in ( |2.3| ) gives 

lt-i{e'n + Jt_i(e-^'^^) = Jt_2(^). (2.8) 

In the UV on the other hand, the current does have a non-trivial monodromy, which 
implies that the right hand side of ( |2.4| ) does not vanish, and eg has a negative real part 
for 1 2; I > 1. Therefore, for \z\ > 1, the fusion relations close with 

Jt_2(e^"z)+Jt_2(e-^-z) = Jt_3+Xt-i, (2.9) 

where now It-i (still with the same proportionality coefficient). In addition, ( ^^ ) 

reads 

j(e*--) _ j(e-*^2) = Ji_i. (2.10) 

In both cases, it is easy to show that these relations are entirely equivalent to those 
of the previous section. 

We thus see that all the monodromies are determined by the fusion relations, together 
with the relation between the current and the partition function ( p.4| ). Knowing these 
relations is like knowing the monodromies, from which the differential equation follows. 
The current is then quickly identified as the solution that is regular at ^ = 0, as a result, 
again, of duality. 

3. The other solutions 

3.1. The other solutions are quasiparticle densities 

The physical significance of the other solutions to the differential equation is easy 
to establish when t is integer. The ek arise in the Bethe ansatz as the energy it takes 
to create a single quasiparticle, and their derivatives correspond to the densities of these 
quasiparticles (or their holes). 

These quasiparticles are familiar from the study of the sine-Gordon model, and consist 
of a soliton and antisoliton, along with t — 2 breathers, which are soliton-antisoliton bound 



12 



states. In the presence of the voltage, there are no breathers in the ground state at zero 
temperature, i.e. only "breather holes". The energy ek required to create a k^^ breather 
of mass mfc = 2 sin ^ 2(t-i) ) computed by the same technique as the one used to 

compute the current. The density of holes is simply related by 2Tip^{6) = ^ oc z^^, 
where, in the notations of 0, the rapidity 9 of the particles is defined by the requirement 
that ek oc m^e^ at large 9. It is related to u via u = e"^"'"^"^, A is the Fermi rapidity; 
z is still related with u by ( |1.5| ), the UV region corresponds to Re{9) > A, the IR to 
Re{9) < A. The densities are precisely related to the Xj. we introduced earlier through the 
normalization 



H it-l)V, 



2np1 = ' ' Ik k = l...t-l. (3.1) 

mt 

Here V , the physical voltage, is proportional to e"^ (see for the exact value of this 
coefficient). 

For /c = t — 1, the meaning of formula ( |3.1|) is that 2txp\_y = 27rp^. In fact, for 9 > A, 
Pa = 0, but for 9 < A, ps = PaT so formula ( |3.1| ) gives this density too. 

For t = 2, the particles do not interact, so there are no breathers. The density of 
states for free particles means that simply 27rpjj; = e^, as is easy to establish for a free 
theory. In terms of the variable z, this goes as which is precisely the second solution 
z~^/'^F (O, ^'i ^'i z) of the hypergeometric equation corresponding to ( p..8| ), as illustrated in 
( p..lO|) . This agrees with the series expansions: from ( |1.8| ), we see that in the region \z\ > 1, 



J(e'''''z) - I(z) = Vn} (-1) « sm — — / - ^ 

^ ^ v; V ; V 2 / r 71 + 1 r 3/2-n/2 V4 



All terms here vanish except the term n = 1, giving X (e^^'^z) — I{z) = —^1727 Siud thus 
1i{z) = — Using that — ^,e^ = ^ for t — 2, this gives rise to 2Tipi = e^, as 
expected. 

The fact that the Xk are related to quasiparticle densities gives a physical significance 
to the radius of convergence for the expansions. In terms of 9, the circle of convergence 
|z| = 1 intersects the real axis for 9 = A. This follows from the fact that A is the Fermi 
rapidity, where some quantities are expected to be singular; for instance, the density for 
solitons vanishes for 9 > A but is non zero for 9 < A. 

Notice that, unlike the current, the other solutions have a singularity at the origin 
as well as at ^ = 1 and z = 00. This has a simple physical interpretation in terms of 
the operators controlling the behavior of these quantities near the IR fixed point. As 
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mentioned in the introduction, the current is controUed only by the tunnehng operator 
cost$ ($ the dual boson); the densities, like thermal properties, depend also on the local 
conserved quantities which are polynomials in the derivatives of $ 

For completeness, we finally would like to give the expansions for the quantities Zbsg 



of [|T0[; one has 



T^lnZBSGiz,±tV/2T) = 

+ V- -S^^jd ^t) e ( -in2it-l)A \ ' 

2t2y^^^ r(n) r[f + n(i-i)] 2zcos^ r V 

3.2. Representing the densities with the curve 

The Ik and hence the densities can also be represented by integrals on the same 
hyperelliptic curve as the current in ( p..4D . This is very similar (and sometimes identical) 
to what happens for the higher-spin magnetizations in the Kondo problem, as discussed 
in 0J^. For simplicity, we will discuss the case t odd only, for which it is convenient to 
make a change of variable x — > -, so the current can be written 

^iu J Co \Jx^-^ + 1 - u^x^ Jco y ' 

Here the contour starts at infinity, loops around the branch point on the positive real axis, 
and goes back to infinity. 

In general, the integrand has t branch points, which, for large u are approximately 
located at Xk = ■u-2/tg2i7rA;/t^ ^ _ 0, ... ,t — 1. One branch point is on the positive real 
axis, and at large u the others are symmetrically distributed in the complex plane. As a 
function of u, the current has singularities when z = 1, an equation with 2{t — 1) solutions. 
As u loops around one of the singularities, z loops around 1, and monodromies are picked 
up. Of course, as u goes around one of the singularities, the different branch points of the 
integrand are exchanged. As a result, or, more simply, by using formula ( |1.17|) , one can 
show the following. The dimensionless reduced densities for breathers with t — k even can 
be written simply as 

Tt-iM = / ^, (3.4) 

•^Cp^t—p y 

where Cp^t-p loops around the singularities Xp and Xt-p, i.e. the (p + lY^ cycle of the 
Riemann surface of the integrand, as shown on the figure { Cq = Ct can be considered the 
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C_23 






c_o 

C_14 



Fig. 1: A choice of cuts and contours in the complex x plane for t = 5. The 
integral around Co gives the current, while integrals around C23 (resp. C14) 
give Xi (resp. X3). 

first cycle). For g = ^ for instance, we thus see that the integral around one cycle is Xi, 
the density of the first breather, while the integral around the second cycle is the current. 

It is also possible to express the dimensionless densities for breathers with t ~ k odd 
by choosing different cycles. For instance, by writing 

one can show that Xt-i [e™z) reads as ( |3.4| ) but with the integral taken along a contour 
Cio that loops around xi and xq. More generally, one has 

=It-2p-i [e^'^z] = [ (3.5) 

as illustrated in figure 2. 

Similarly, Jt-2p-i {e~'^'^z) reads as (|3.5| ), but with a complex conjugate contour. As 
a result, fusion relations take a very simple form. For instance, one has Its {e^^z^ + 
Xi_3 (^e~^'^z) = Xt-2{z) +2t-4:{z), is illustrated on figure 3. 

Also, the truncation of the fusion relations 2t-i (e^^z) +2t-i {e~^'^z) = 2t-2{z) can 
be illustrated graphically in a similar fashion. 

To conclude, we illustrate these considerations by the case t = 2. We now have only 
two branch points, and the integral of the differential form along the unique cycle is (going 
back to the original x variable) 

dx __TT 

4 J^^ uVx + a;2 - ^2 2u ' 
15 



CJ2 




c: 



C 34 



Fig. 2: An other choice of cuts and contours in the complex x plane for 
t = 5. The integral around Cq gives the current, while integrals around C12 
(resp. C34) give X2 (resp. X2). 




Fig 



2ie '""z 



3: Graphical representation of the fusion relation X2 ( 

) =Isiz)+Ii{z) fort = 5. 
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where xq and xi are the two branch points. On the other hand, 
SO Xi (e*^^) = — 7r/(2M) as required. 



4. Conclusion and Questions 

This paper we hope is a first step in finding new ways to solving integrable models. 
The case we have discussed can be thought of as the deformed SU (2)i case (when t = 1 the 
bulk conformal field theory is the SU {2)i WZW model). This has been well understood 
by many approaches, and so our results do shed light on the picture but don't result in 
the computation of any new quantities. In the SU{2)k case (the multi-channel quantum 
wire and multi-channel Kondo problem), the zero-temperature curves are known and 
the multi-channel Kondo problem fairly well understood (see and references within). 



The functional relations can easily be derived for most of the quantities, but the crucial 
one involving the current (the analog of (|2.4| )) is not yet known. In the SU{N) case,, the 
problem is even murkier. The SU{N) Kondo problem is integrable [jl^, but only in a 
special case has the zero-temperature magnetization been calculated. Moreover, it is not 
even known which if any SU{N) generalizations of the boundary sine-Gordon problem are 
integrable. We hope that our approach will be useful in understanding these problems. 

More formally, observe that, although we were motivated by the study of the boundary 
sine-Gordon model, all the results obtained do have a well defined meaning within the c = 1 
conformal field theory itself, and its integrable structure. It is tempting to speculate that 
monodromy and duality should play a crucial role in the study of this integrable structure, 
in particular that there should be hyperelliptic curves "naturally" associated with rational 
conformal field theories. More progress in understanding the general picture is clearly 
needed here. 

Finally, we find it intriguing that the curves are so simple and continuous in t, given 
that the differential equations change dramatically as t is varied. It is hard not to believe 
that there is a simple and deep reason underlying this, but we do not know what this is. 

This research was supported by the NSF (DMR-9802813), the DOE (Outstanding 
Junior Investigator award) and by the Sloan Foundation (P.F.); and by the NYI program 
(NSF-PHY-9357207) and DOE grant DE-FG03-84ER40168 (H.S.). 
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